Enhanced the thermal Entanglement in Anisotropy Heisenberg XYZ Chain 
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The thermal entanglement in Heisenberg XY Z chain is investigated in the presence of external 
magnetic field B. In the two-qubit system, the critical magnetic field Be is increased because of 
introducing the interaction of the z-component of two neighboring spins Jz- This interaction not 
only improves the critical temperature Tc, but also enhances the entanglement for particular fixed 
B. We also analyze the pairwise entanglement between nearest neighbors in three qubits. The 
pairwise entanglement, for a fixed T, can be strong by controlling B and Jz. 
PACS: 03.65. Ud, 03.67. -a, 75.10. Jm 

rn" 

O ■ I. INTRODUCTION 

o 

(N 



^3 



Entanglement is an important resource in quantum information [1] . The ideal case which Quantum computing and 
quantum communication are put into use is to find entanglement resource in solid system at a finite temperature. 
I The Heisenberg model is a simple but realistic and extensively studied solid state system [2] [3] . Recently, Heisenberg 
] interaction is not localized in spin system. It can be realized in quantum dots [4], nuclear spins [5], cavity QED [6] 
' [7]. This effect Hamiltonian can be used for quantum computation [8] and controlled NOT gate [7]. The thermal 
^ I entanglement in isotropic Heisenberg spin chain has been studied in the absence [9,10,15] and in the presence of 
. an external magnetic B [9,10,14]. The entanglement of two-qubit isotropic Heisenberg system decreases with the 
increasing T and vanishes beyond a critical value [9,10], which is independent of B. Pairwise entanglement in 
, 7V-qubit isotropic Heisenberg system in certain degree can be increased by increasing the temperature or the external 
field B [9]. Anisotropic Heisenberg spin chain has been investigated in the case of B = [10] and B ^ [11]. For 
, a two-qubit anisotropic Heisenberg XY chain, one is able to produce entanglement for finite T by adjusting the 
' magnetic field strength [11]. However, the entanglement by increasing T or B, in two-qubit anisotropic Heisenberg 
XY chain [11] or in 7V-qubit isotropic Heisenberg chain [9], is very weak. How to produce strong entanglement is 
worthy to study. 

On the other hand, we have not find the work about two-qubit or the iV-qubit anisotropic XYZ Heisenberg chain 
in the presence of magnetic field. Although the 7V-qubit Heisenberg chain has been studied [12] [9], in Ref. [12] the 
authors studied the maximum possible nearest neighbor entanglement for ground state in a ring of iV qubits, and in 
^ [9] they just investigated the case of isotropy N qubits Heisenberg chain. In this paper, we study the entanglement 
O^' of two-qubit anisotropic Heisenberg XYZ chain and the pairwise entanglement of three-qubit anisotropic Heisenberg 
^ I XYZ chain. Introducing the interaction of the z-component of two neighboring spins not only improve the critical 
temperature Tc but also enhance the entanglement for fixed B and T in particular regions. In the case of anisotropic 
three-qubit Heisenberg XYZ chain, the effect of partial anisotropy 7 make the revival phenomenon more apparent 
than in two-qubit chain; for a fixed T, one can obtain a robust entanglement by controlling B and Jz- 
The Hamiltonian of A^-qubit anisotropic Heisenberg XYZ model in an external magnetic field B is [11] 
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- ^[J.af af+i + J,af f7,:Vi + + B{a! + af+i)], (1) 

i=l 

where 'o^j ~ {a^ , a^- , a^^) is the vector of Pauli matrices and Ji{i ~ x,y,z) is real coupling coefficient. The coupling 
coefficient Ji of arbitrary nearest neighbor two qubits is equal in value. For the spin interaction, the chain is said to 
be antiferromagnetic for > and ferromagnetic for < 0. 

For a system in equilibrium at temperature T , the density operator is p = Z^^ c'xjp{—H/kBT), where Z = 
Tr[exp(— iJ/fcsT)] is the partition function and ks is Boltzmann's constant. For simplicity we write ks ~ 1. Entan- 
glement of two qubits can be measured by concurrence C which is written as C = max(0, 2max{Ai} — X^iLi ^i) [16] 
[17] [13], where Xi is the square roots of the eigenvalues of the matrix R — pSp*S, p is the density matrix, S = af^^al^ 
and * stand for complex conjugate. The concurrence is available no matter what p is pure or mixed. 
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II. TWO-QUBIT HEISENBERG XYZ CHAIN 



Now, we consider the Hamiltonian for anisotropic two-qubit Heisenberg XYZ chain in an external magnetic field 
B. The Hamiltonian can be expressed as 

J B 

H = J{a+a2 + a^a+) + Jj{a+a+ + a^a^) + y<T?(T^ + + (^2) (2) 

where = \{'^'^ ^ iay) is raising and lowering operator respectively, and J = ^^ii^ ^ = TJ^e parameter 7 

(0 < 7 < 1) measure the anisotropy (partial anisotropy) in XY plane. When the Hamiltonian of the system has the 
form of Eq.(2), in the standard basis {|00), |01), |10), |11)}, the density matrix of the system can be written as 

(Ml u \ 
u; z I , . 

Q z w Q ]■ 
V M2 / 

These nonzero matrix element can be calculated through 

Mi=Tr(|00)(00|p),M2-rr(|ll)(ll|p), 

w = Tr(|01)(01|/9), w = Tr(|00)(ll|p), z = Tr(|01)(10|p). (4) 

The square roots of the eigenvalues of the matrix i? are Ai ^2 = jwit zj, A3_4 = |^uiM2±w|. Therefore, we can calculate 
the concurrence. 

The eigenvalues and cigcnstates of H are easily obtained as = (— ^ ± J)|^'^), = ± r7)|E^),with 

the cigcnstates [f^) = ^(|01) ± |10)), = --^=i_[(?7 =f -B)|00) ± J7III)], where 77 = y^W+JlW- One can 

notice that the eigenstates are the same as the case of = [11]. Because the basises |01) and |10) are the two 
degenerate eigenstates of (jI(J2 with eigenvalue —1, hence the superposition of the two degenerate states |01) and |10) 
still is the eigenstate of crfcrl, that is, |\E'^) is the eigenstate of = as well as that of 7^ 0. The same reason 
account for |S='=) both as an eigenstate of Eq.(2) and as that of the case of Jz = . From Eq.(4), tracing on the 
eigenstates, we obtain the square roots of the eigenvalues of the matrix R 



Ai,2 = Z ^ 
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A34 = Z-^e-^\Jl + {^smhf3r]Y T — sinh/???!, (5) 
V ^ ^ 

where the partition function Z = 2(e^2T cosh/Jry + cosh /3 J). Because the concurrence is invariant under the 
substitutions J ^ —J and 7 — > —7 [11], we will consider the case J > and < 7 < 1. But with substitution 
Jz — > —Jz the the concurrence is variant. We choose Jz > 0, and we will state the reason later. 

We first review the circumstance of anisotropic Heisenberg XY chain, which is analyzed in [11]. At T = 0, exist a 
critical magnetic field Be- As B cross Be, the concurrence C drops suddenly then undergoes a "revival" for sufficient 
large 7. However, we noticed that Be decrease with the increasing of the anisotropic parameter 7. Although with 7 
increasing the critical temperature Tc is improved, the entanglement, when temperature is in the revival region, is 
very weak. 

With 7 = 0.3, we show the concurrence as a function of B and T for two values of Jz in Fig. 1. For Jz = (Fig. la) 
corresponding to the circumstance of anisotropic Heisenberg XY chain [11], one can observe a revival phenomenon 
and the weak entanglement in revival region. For the convenience of representation, we define the main region in 
which concurrence C keeping its constant and maximal value. Comparing Fig. 1 (a) with (b), we find that with the 
increasing of Jz, the main region is extended in terms of B and T, i.e., the critical magnetic field Be is broadened 
and the critical temperature Tc in main region is improved. That is to say, the range of concurrence C keeping its 
constant and maximal is extended in terms of B and T, so we can obtain strong entanglement in the extended range. 

We can understand the efi^ect of Jz on Be from the case of T = 0. For T = under the condition of Jz < J, C can 
be written analytically as 

ioi 1] < J + Jz 

C{T = 0) = { (1 - Jl/v)/2 iovj] = J + Jz (6) 
/7/r/ ior T] > J + Jz 
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The parameters J, rj and 7 arc independent of in the case of two interacting qubits. Comparing Eq.(6) with Eq.(6) 
of Ref. [11], we can see clearly that if Jz is positive, makes the intersection points of piecewise function shift. In 
this paper, we consider the case oi Jz > 0. Fig. 2 shows the concurrence at T = for three vahies of positive Jz- 
It show clearly that concurrence drops sharply at a finite value of magnetic field B, which is called critical magnetic 
field i?c, at which the quantum phase transition occurs [11]. But with the increasing of Jz, is increased. The 
interaction of the z-component of two neighboring spins Jz causes a shift in the locations of the phase transitions. 
Namely, the presence of positive Jz increases the region over which the concurrence C attains its maximum value. This 
result means that in larger region of B and T we can obtain stronger entanglement. The effect of Jz is different with 
that of 7 on changing Be. In the case of = [11], although with the increasing of 7 the critical temperature T,, 
is increased, the larger the values of 7, the smaller the critical magnetic field Be- Here, introducing the z-component 
interaction of two neighboring spins not only extends critical magnetic field Be but also improves critical temperature 
Tc and the entanglement (we will further show it in Fig. 3). 

Let us consider concurrence changing with temperature for different values of Jz in a fixed B (B = 1.1). We plot 
it in Fig. 3 with 7 = 0.3. We notice that existing a critical temperature Tc at which the entanglement vanishes. 
Obviously, Tc is improved monotonously with increasing of Jz- Under the condition = (corresponding to XY 
model [11]), the concurrence exhibit a revival phenomenon, but the maximal values of entanglement in both area are 
small. If introducing the Jz, the critical external magnetic field Be become larger so that B = 1.1 is less than Be 
(the critical magnetic when Jz = 0.2,0.5 or Jz = 0.9), thus we observe the maximal value of entanglement 1. In 
the temperature range < T < 1.725 , the larger the stronger entanglement. Therefore, Jz not only improve the 
critical temperature Tc, but also enhance the entanglement for particular fixed B and 7. 

III. THE PAIRWISE ENTANGLEMENT IN THREE QUBITS 

The calculation of pairwise entanglement in N qubits is very complicated due to the anisotropy in Heisenberg XYZ 

chain. Here we just calciilate the pairwise entanglement in three qubits to show the effects of Jz . We now solve the 
eigenvalue problems of the three-qubit XYZ Hamiltonian. We list the eigenvalues and the corresponding eigenvectors 
as follow 

£^1,2 = -J-:^+B: = ±i(l T ^)|110) + ^1101) T i(l ± ;^)|011)), 

T 1 2 

^3,4 = J + f - S ± r?_ : |$3,4) = ,^ ^ -2B- J± r?_)|000) + J7 E Tf"|110)]; 

2 V 2»?- [V- ± {Jz -2B- J)] ^0 

E,,e = -J-^-B: J^s.e) = ±^(1 T ^)|010) + -^JlOO) T ^(1 ± ;^)|001)); 

T 1 2 

Er,s = J+f+B±n+: J^r.s) = , . ^.^ ^A {J. +2B-J± r?+)llll) + J7 ^ T"1010)]. (7) 

2 ^2r/+ [r]+ ± {Jz +2B- J)] 



where t]± = \J{Jz — J ^ 25)2 _|_ 3( j^)2^ y jg ^-j^g cychc right shift operator [15]. The reduced density matrix of two 
nearest- neighbor qubits in N qubits system also has the form of Eq.(3). Employing Eq.(4) and tracing on the basis 
of eigenstates shown in Eq. (7), one can get the density matrix [i-^, z,v, then further obtain the concurrence. 

Here we do not write the expressions of Aj because it is very long. We will directly plot some curves to show the effect 
of Jz on enhancing entanglement. 

Fig. 4 show concurrence as a function of B and T with 7 = 0.3, Jz = 0.9 and J = 1.0 in three-qubit XYZ 
Heisenberg chain. We see that with the same 7 = 0.3, the effect of partial anisotropy 7 make the revival phenomenon 
more apparent than in two-qubit chain. When B = A in Fig. 1, the largest critical temperature Tc produced by 7 
is about 1.0 (Fig. la); due to the restrain of Jz the maximum temperature only caused by 7 is below 0.8(Fig.lb). 
However, in three-qubit system if B = 4 with the same set of parameters, comparing Fig. lb with Fig. 4, the critical 
temperature Tc in revival region almost equal to 1.8. The stronger effect of 7 implies that if we aim to obtain strong 
entanglement we can decrease 7 properly and increase Jz, otherwise increasing 7 can make the revival phenomenon 
more evident. Of course, the coupling constant Jz also increase magnetic field Be and expend the region of concurrence 
keeping constant in terms of B and T as it do in two-qubit (for the limited of the page,we do not plotted here) . 

For T = 0.6, Fig.5 show concurrence as function of B and Jz- There is no entanglement for B = 0, which corresponds 
with Fig. 4. If Jz is below a certain value, in case of Fig. 5 the value is about 0.2, the entanglement appears in one 

area corresponding to the "revival" [11] on condition that the magnetic field is larger than a certain value, and the 
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certain value of B is increased with the enhanced of J^. But, if is larger than 0.2, there are two areas appearing 
entanglement, and the entanglement appearing in the lower range of B can be much stronger than that in higher 
magnetic field. In the lower range of B, for a certain B, the large the large concurrence. Thus, in the A''-qubits 
XYZ system, for a fixed T, one can obtain a robust entanglement by controlling B and J^. 



IV. CONCLUSION 



The thermal entanglement in anisotropic XYZ Heisenberg (;hain is investigated. Through analyzing the two-qubit 
system, we find that with the increasing of Jz, the critical magnetic field Be is increased; the coupling along Z not 
only improves the critical temperature Tc, but also enhances the entanglement for certain fixed B. We also analyze 
the entanglement between two nearest neighbors in three qubits and find that the effect of partial anisotropy is more 
evident than it do in two-qubit system. The pairwise entanglement exhibit a interesting phenomenon. For certain 
fixed B , if the coupling constant Jz is small, the pairwise entanglement only exists in relative strong magnetic field 
B and the entanglement is weak. By increasing Jz, in lower range of B , one can obtain a strong entanglement. 
Therefore, interaction constant of the z-component of two neighboring spins Jz play important role in enhancing 
entanglement and in improving the critical temperature. 

This work was supported by Ministry of Science and Technology of China under Grant No.2100CCA00700 

The captions of the figure: 

Fig. 1 Concurrence in two-qubit Heisenberg XYZ chain is plotted vs T and B, where (a): Jz = 0, (b): Jz = 0.9. 
For aU plotted J = 1.0, 7 = 0.3. 

Fig. 2 Concurrence in two-qubit Heisenberg XYZ chain vs B at zero temperature for various values of Jz with 
7 = 0.3 and J = 1.0. Prom left to right Jz equal to 0, 0.5, 0.9, respectively. 

Fig. 3 Concurrence in two qubits Heisenberg XYZ chain is plotted vs T . For all plotted J = 1.0, B = 1.1,7 = 
0.3. From top to bottom Jz equal to 0.9, 0.5, 0.2, 0, respectively. 

Fig. 4 Pairwise entanglement in three-qubit Heisenberg XYZ chain is plotted as a function of T and B, where 
7 = 0.3, J= 1.0, Jz = 0.9. 

Fig. 5 Pairwise entanglement is plotted as a function of B and Jz, where T = 0.6, J = 1.0, 7 = 0.3 . 
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